Bloch Structures in a Rotating Bose-Einstein Condensate 
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A rotating Bose-Einstein condensate is shown to exhibit a Bloch band structure without the need 
of periodic potential. Vortices enter the condensate by a mechanism similar to the Bragg reflection, 
if the frequency of a rotating drive or the strength of interaction is adiabatically changed. A localized 
state analogous to a gap soliton in a periodic system is predicted near the edge of the Brillouin zone. 



The Bloch band is crucial to our understanding of the 
behaviors of periodic systems, such as electronic states 
in solids Q, atom dynamics in optical lattices Q, and 
Cooper-pair tunneling in small Josephson junctions 
The underlying physics common to all of these is the 
Bragg reflection that occurs at the edges of the Brillouin 
zone. The resultant Bloch band leads to various interest- 
ing phenomena, such as Bloch oscillations [l| (a particle 
in a periodic potential driven by a weak constant force 
cannot be accelerated indefinitely but oscillates in real 
space) and formation of gap solitons (localized wave 
packets arising from the balance between negative-mass 
dispersion and repulsive interaction); both phenomena 
have recently been observed with a Bose-Einstein con- 
densate (BEC) in an optical lattice 0, El • 

In this Letter, we show that yet another system — a 
BEC confined in a rotating harmonic-plus-quartic poten- 
tial — exhibits a Bloch band structure, and investigate 
the associated novel phenomena. Seemingly this system 
has no periodic structure, but may be considered to be 
a quasi-lD periodic system in the following sense. When 
the rotating frequency of the potential is high, the quartic 
potential oc r 4 , together with the centrifugal one oc — r 2 , 
produces a Mexican-hat shaped potential 01 whose min- 
ima form a quasi- ID toroidal geometry. Under such cir- 
cumstances any perturbation V(9) that is needed to drive 
the system into rotation by breaking the axisymmetry 
of the potential serves as a "periodic" potential since 
V{9) = V(9 + 2tt), where 9 denotes the azimuthal an- 
gle. 

In the present system, we will show that the Bragg 
reflection causes vortex nucleation. Since this process 
occurs adiabatically, in contrast to a method invoking 
dynamical instabilities 0, 0] , we can obtain the vortex 
states without heating the atomic cloud. A stirring tech- 
nique to produce vortices without heating has also been 
proposed in Refs. [ToL Hl| . Another scheme that adia- 
batically nucleates vortices is based on phase engineering 
techniques fl2l H?3 . Ilil ] . We will also show that a localized 
state is generated even with repulsive interactions, which 
is analogous to gap solitons in periodic systems 

We begin by discussing a BEC in a ID ring to under- 
stand the essence of the phenomena. We assume that a 
rotating potential takes the form V(9, t) = ecosn(9— tit). 



In a frame rotating at frequency Q, the potential V be- 
comes time-independent, and the Gross-Pitaevskii (GP) 
equation is given by 
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where energy and time are measured in units of 
h 2 1 (2mR 2 ) and 2mR 2 /h with m and R being the atomic 
mass and the radius of the ring. The wave function is 
normalized as \ipii)\ 2 d9 = 1, and 7 characterizes the 
strength of interaction. 

First, let us consider the noninteracting case (7 = 0). 
Introducing a new variable cf> = e~ ine / 2 ipiD, we transform 
Eq. into the Mathieu equation (—d@ + ecosn9)tfi — 
Etfi, where E = fi + fi 2 /4 plays the role of the "total 
energy" in our Bloch band picture. Then the solution 
satisfies q)(9 + O ) = e^ p - n /^ @0 4>(9), with 6 = 27r/n 
being the period of the potential V , and p an integer. It 
follows that p — VI / 2 may be regarded as the quasimo- 
mentum, indicating that we can move in the quasimo- 
mentum space by changing Q. Since the unit "reciprocal 
lattice vector" is 2n/Qo = n, the Bragg reflection occurs 
only between the states whose angular momenta differ 
by n, 2n, ... (Bragg's law). Consequently there are n 
independent branches corresponding to p = 0, • • ■ , n — 1 
modulo n. Figure ^ shows E and the angular momen- 
tum (L) = —i J 2 ^ ipiY)dgipiT)d9 as functions of Q for the 
first and second Bloch bands with n = 2 and e = 0.1. 
Two independent branches exist in this case, and are 
distinguished by the solid and dotted curves. The Bloch 
band structure in Fig. ^ indicates that if we adiabati- 
cally increase or decrease fi across = 2, the nonvor- 
tex state -0id — l/V2n transforms to the vortex state 
■0id — e 2ie /^/2tt, and vice versa. The time scale of the 
adiabatic change must be much longer than the inverse 
energy gap ~ 

When |7| > e, the system exhibits hysteretic behavior, 
and a loop structure emerges in the Bloch bands as 
shown in Fig. We see that the first (second) band 
forms a loop for 7 = 1(— 1); then the adiabatic nu- 
cleation of vortices is not possible with increasing (de- 
creasing) fi because the nonlinear Landau-Zener transi- 
tion [la] occurs no matter how slow the change of £1. It 
can be shown that the states on the loop have dynam- 
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FIG. 1: (a) E = fi + f7 2 /4 and (b) (L) = -i J ip$ D d e ipi O d0 
of the first (red curves) and the second (blue curves) Bloch 
bands for the noninteracting (main panels) and interacting 
(right panels) BECs in a ID ring described by Eq. Q with 
V{9) = 0.1 cos 26. The solid and dotted curves refer to two 
independent branches. 

ical instability that breaks the symmetry of the system 
|^id(#)| = |V'id(^ + ©o)|- This instability is similar to 
the one that forms gap solitons in a BEC in an optical 
lattice [lfij |. 

We now consider a BEC in a harmonic-plus-quartic 
potential with a rotating stirrer. We assume that the 
axial trapping energy huj z is much larger than the other 
characteristic energies (a tight pancake-shaped trap) so 
that the system is effectively 2D. The external poten- 
tial in the corotating frame of reference takes the form 
V(r) = r 2 /2 + Kr 4 /4: + er 2 cos29, where if is a con- 
stant and the third term e(x 2 ~ y 2 ) is a stirring potential 1 \^\ 4 /2}d6 = m 2 -Q.m + 7/(4tt) for tp^> = e im6 /V2^. 
which can be produced by laser beams propagating in the 
z direction. Such an anharmo nic p otential was theoreti- 
cally considered in Refs. [71 l l7L and has recently been 
realized experimentally 19]. We normalize the length, 
time, energy, and wave function by do = (h/muij_), u)J_ , 
fiu>±, and VJf/do, respectively, with uo±_ and N being 
the frequency of the radial harmonic trap and the num- 
ber of atoms. The wave function is then normalized as 
J \ fp2B\ 2 dr — 1. The time-dependent GP equation in the 
frame corotating with the stirring potential at frequency 
£1 is given by 



fective strength of interaction in 2D |2(| with a being the 
s-wave scattering length. 

When ft>l, the system described by Eq. J5J) can be 
approximated by a quasi- ID ring |7j. Assuming ip2T> — 
f(r)e iee /V2~jr, we find that an effective potential for the 
radial wave function is given by £ 2 / (2r 2 ) + r 2 / 2 + Kr 4 / 4. 
This potential has a minimum at r ~ (e/K) 1 ^ for 
i 1, and the effective frequency around it is lu c s ~ 
^(iK) 1 / 3 . If uj c g is much larger than other character- 
istic frequencies, the dynamics of f(r) can be ignored. 
After normalization of the time by p~ 2 = J r\f\ 2 r~ 2 dr, 
the equation of motion for 9 reduces to Eq. (Q, where 7 
is given by 2gp 2 J r\f\ 4 dr. Thus, the system described by 
Eq. J3J) exhibits the quasi- ID circular flow for large an- 
gular momentum 0, and we expect that a Bloch band 
structure emerges. We show below that this is true even 
for fi ~ 1 and I ~ 1. 

Figure [21 (a) shows the angular momentum (L) = 
—i J , 02D^V'2Drfr and the density and phase pro- 
files (insets) of the noninteracting stationary states of 
Eq. (0) 21]. We find that the behavior of (L) is similar to 
that in Fig. (b) . When we start from the ground state 
with (L) — (1-1), two vortices enter at fi ~ 1.57 (1-2), 
producing a doubly-quantized vortex (1-3) [l7|. Increas- 
ing the rotation frequency £1 further, we can nucleate two 
more vortices in the condensate [(1-4) and (1-5)]. If we 
follow the branch starting from (L) = 2 (2-1), the two 
vortices escape out of the condensate with increasing f2 
[(2-2) and (2-3)], and then four vortices enter at fi ~ 1.7 
(2-4). It should be noted that the Bloch band picture 
holds even for the nonrotating gas, while the ring-shaped 
profile manifests itself only if the gas is rotating [cf. (1-1) 
and (2-3)]. 

The frequencies at which the Bragg reflection oc- 
curs in Fig. |3 (a) are different from those in Fig. ^ 
The difference arises from the dispersion relation of the 
vortex states in 2D, i.e., the relation between the en- 
ergy and angular momentum. In the ID case, the en- 
ergy is given by E™ = J{^*(-d% + tftdg)^ 



2D 



dt 



l\7 2 + l n-^- + V(r)+g\i: 2D \ 2 )i, 2D7 (2) 



where g = [uj z /(2'KUJj_)] 1 ^ 2 4:TTNa/do characterizes an ef- 



The Bragg reflection between -0^P and 4'm+2 occurs 
when = E™ +2 , i.e., Q = 2m + 2. In the 2D 

case, we employ a variational wave function ip 2 ^ = 
(7r|m|!)- 1 / 2 d~ l '" h Vl m l exp[-r 2 /(2c4) + im9] to obtain 
E™ = (\ m \ + l)(d m 2 + d 2 m )/2 + K(\ m \ + l)(\ m \ + 2)d 4 m /4- 
Vlm + g(2|m|)!/[2 2 l m l +2 7r|m|! 2 d 2 „], where the variational 
parameter d m is determined from dE 2 ^ / dd m — 0. When 
K = 1 and g = 0, the condition for Bragg reflections to 
occur between ^ D and ^ D +2 {E™ = E™ +2 ) for m = 0, 
2, and 4 is satisfied at f2 = 1.58, 1.83, and 2.03, respec- 
tively, which are in good agreement with the numerical 
results in Fig. El (a) (corresponding to the points at which 
the red and blue curves cross with (L) changing by ±2). 

The energy gap can also be obtained by a variational 
method. The matrix element of the stirring potential, 
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FIG. 2: (a) Angular momentum (L) versus rotation frequency 
SI of the stationary states of the GP equation ||2J with g — 0, 
K = 1, and e = 0.02. The images (l-l)-(l-5) [(2-l)-(2-5)] cor- 
respond to the branch indicated by the red (blue) curve, (b) 
Time evolutions with A = 0.03, where SI is linearly changed 
from 1.52 to 1.6 (red curve) and from 1.62 to 1.54 (blue curve) 
during < t < 600. The initial state is the ground state, and 
e is linearly ramped up from to 0.04 during < t < 200, kept 
at 0.04 during 200 < t < 400, and ramped down from 0.04 to 
during 400 < t < 600 to suppress nonadiabatic disturbances. 
The size of the images is 5 x 5 in units of (h/muj^) 1 ^ 2 . 



for example, between i/jq and ■02 * s calculated to be 
1/ ip% D *ipo B sr 2 cos29dr\ ~ 0.46e for K = 1. The energy 
gap at SI ~ 1.57 is then given by ~ 0.92e, which also 
agrees very well with the numerical result ~ 0.9 le. 

The Bloch band structure in Fig. 0(a) indicates that if 
we prepare the nonvortex state at SI < 1.55, and adiabat- 
ically increase SI to above SI ~ 1.6, we obtain a doubly- 
quantized vortex state by following the red curve. Inter- 
estingly, we can nucleate vortices also by decreasing Q, 
that is, by preparing a nonvortex state for 1.6 < SI < 1-7 
and decreasing SI to below SI ~ 1.55 by following the 
second Bloch band (blue curve) . Figure (b) illustrates 
these adiabatic processes. In order to examine the ef- 
fect of dissipation on the adiabatic processes, we replace 
i with i + A on the left-hand side of Eq. J5J), where a 
constant A is taken to be 0.03 j^- We prepare the non- 
vortex ground state with g = 0, and gradually ramp up 

from 1.52 to 1.6 [red curve in Fig. 13(b)], or ramp it 
down from 1.62 to 1.54 (blue curve). As expected, two 
density holes with phase singularities come from infinity, 



FIG. 3: (a) The g dependence of the first and second bands 
(left and right panels) with K — 1 and e = 0.02, where the 
curves are offset by ±g for clarity. The dynamically unstable 
regions are indicated by the red curves. The inset shows a 
localized stationary state - "gap soliton"- formed near the first 
band at g = 0.8 and f2 = 1.55. (b) Time evolutions with 
A = 0.03 and K = 1, where g is linearly ramped from —1.2 to 
0.5 with Q — 1.58 [red curve, corresponding to the red arrow 
in (a)], from 1.5 to —0.3 with SI = 1.56 (blue curve, blue 
arrow), and from to 1 with £1 = 1.55 (green curve, green 
arrow). A small perturbation to break the axisymmetry is 
added to the initial state. The change of e in the red and 
blue curves is the same as in Fig.|2](b) , and e = 0.04 is kept 
during 200 < t < 600 in the green curve. 



and unite to form a doubly-quantized vortex. We note 
that the behaviors in Fig. (b) are almost the same as 
for the dissipation- free case (A = 0), indicating that the 
process is robust against dissipation. The energy gap 
between the first and second Bloch bands at ~ 1.57 is 
~ 0.036 for e = 0.04. According to the Landau-Zener for- 
mula, the transition probability between the energy gap, 
i.e., the probability that the nonvortex state remains so, 
is below 1% in the situation in Fig.[5](b), in agreement 
with our numerical result. 

In the case of g ^ 0, the interaction bends the Bloch 
bands as shown in Fig. [21 (a), yielding hysteresis as in the 
ID case. We note that the region of SI in which the Bragg 
reflection occurs is shifted due to the interaction. The 
shift is attributed to the difference in interaction energies 
between the nonvortex and vortex states. To show this 
we again employ the variational wave function tp^. We 
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find that the condition for the Bragg reflection to occur 



E, 



2D 



, is satisfied at 



between tpQ and ip. 
Q = 1.53 for g — 1.6 and SI = 1.64 for g = —1.6, in 
qualitative agreement with the shifts shown in Fig. [3] (a). 
Since the interaction term of is independent of to, 
the shift does not occur in ID (see Fig.^J. 

The shift in the position of the Bragg reflection im- 
plies the possibility of the adiabatic nucleation of vor- 
tices by changing the st reng th of interaction using, e.g., 
the Feshbach resonance ;23j, at fixed f2. For example, at 
il ~ 1.58, the first Bloch band has angular momentum 
(L) ~ for g < — 1 [left panel in Fig. 01(a)], which con- 
tinuously increases to (L) ~ 2 with an increase in the 
interaction to repulsive (red arrow). Similarly, there is a 
region in which (L) changes from to 2 as g decreases 
from positive to negative [blue arrow in the right panel 
in Fig. (a)] . No dynamical instability is present on the 
red and blue arrows. Figure|3|(b) illustrates these results, 
where the initial state is the nonvortex ground state, and 
the strength of interaction is gradually changed with il 
held fixed (red and blue curves). We find that the two 
vortices are nucleated in a manner similar to that in Fig. [21 
(b) . The repulsive-to-attractive case is particularly inter- 
esting because the attractive interaction is usually con- 
sidered to hinder vortex nucleation [lsj . 

The adiabatic theorem breaks down when the path en- 
ters the region of dynamical instability. A notable exam- 
ple is shown as the green arrow in Fig.|3](a). For g > 0.6 
at SI ~ 1.55, a twofold symmetric state [like the inset 
(1-2) of Fig. [3 (a)] becomes dynamically unstable against 
localization into the state shown in the inset of Fig. (a) . 
This localization is due to the interplay between repulsive 
interaction and negative-mass dispersion around the edge 
of the Brillouin zone. The localized state can therefore 
be regarded as an analog of the gap soliton in a periodic 
system 0. The green curve in Fig. El (b) demonstrates 
the time evolution corresponding to the path shown as 
the green arrow in Fig. [31 (a) ■ The twofold symmetry as 
shown by the inset at t = 300 is broken in the course of 
vortex nucleation, giving rise to two localized states as 
shown by the insets at t = 400 and t = 560. 

Figure |2| (a) shows that a doubly-quantized vortex 
state with (L) ~ 2 is dynamically unstable (stable) in 
most of the region with g < (g > 0). This instability 
also gives rise to a localized state that undergoes center- 
of-mass rotation. We find from the numerical analysis 
that a doubly-quantized vortex state with g < evolves 
into the localized state, followed by the oscillation be- 
tween these two states in a manner similar to the split- 
merge oscillations reported in Ref. [24|. 

In conclusion, we have shown that the Bloch band 
structure arises in a BEC confined in a harmonic-plus- 
quartic potential with a rotating drive, which enables us 
to nucleate vortices adiabatically. The physical mecha- 
nism of the vortex nucleation is very similar to the Bragg 



reflection at the edge of the Brillouin zone. Interestingly, 
we can nucleate vortices not only by increasing the stir- 
ring frequency but also by decreasing it, or by changing 
the strength of interaction at a fixed stirring frequency, 
which may be called a "Feshbach-induced Bragg reflec- 
tion". The adiabatic processes are robust against dissi- 
pation due, e.g., to a thermal cloud. Spontaneous local- 
ization of the rotating cloud is predicted, which signals 
the existence of a gap soliton. 
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